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Abstract 
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I. INTRODUCTION 



It is observationally supported that the current expansion of the universe is accelerat- 
ing l|, l2|. The scenarios to account for the current accelerated expansion of the universe fall 

111 ]. One is to introduce "dark energy" in 



into two broad categories 



0.y.B.y.B.y.y.Q 



the framework of general relativity. The other is to study a modified gravitational theory, 



e.g., F{R) gravity, in which the action is represented by an arbitrary function F{R 



scalar curvature R (for reviews, see 



0,0, Is 



10 



Ih and for a new approach, see 



12|). 



the 



131] imply that the ratio of the effective 



On the other hand, various observational data 
pressure to the effective energy density of the universe, i.e., the effective equation of state 
(EoS) Wcs = Poff /Peffi may evolve from larger than —1 (non-phantom phase) to less than — 1 



(phantom phase IJ, ll5|). Namely, it crosses —1 (the phantom divide) at the present time 
or in the near past. A number of models to realize the crossing of the phantom divide have 
been proposed (for a detailed review, see {4]). 

There are also several studies for the crossing of the phantom divide in the framework 

of F{R) gravity , Q, Q, Q, Q, y, I21I . I22P (for related works, see |25|). An exp licit 

model with realizing a crossing of the phantom divide has been constructed in Ref. 26|] and 



its thermodynamics has been examined 27|]. Moreover, in Ref. 



28 



it has been illustrated 



that multiple crossings of the phantom divide can occur in F{R) gravity as the scalar field 
theories such as an oscillating quintom model {29 or^a quintom with two scalar fields [s"" 
in the framework of general relativity (see also 31 



32|). 



In this paper, we study a model of F{R) gravity in which a crossing of the phantom divide 
can be realized by taking into account the presence of cold dark matter. We demonstrate 
the behavior of F{R) gravity around a crossing of the phantom divide. In our previous 



work 



261 ]. an analytic solution of F{R) gravity to realize a crossing of the phantom divide 



without matter has been derived. In this work, as a further investigation, we examine a 
solution of F{R) gravity to realize a crossing of the phantom divide with cold dark matter. 
We use units oi = c = h = 1 and denote the gravitational constant SnG by = Stt/Mpi^ 
with the Planck mass of Mpi = G^^/^ ^ 12 x lO^^GeV. 

The paper is organized as follows. In Sec. II, we explain the reconstruction method of 



^ The equivalence between F{R) gravity and the scalar-tensor theory has been indicated in Ref. [23|. The 
crossing of the phantom divide in scalar-tensor theories has been investigated in Ref. 24 1. 
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F{R) gravity proposed in Ref. 33j. In Sec. Ill, we study a model of F{R) with realizing 
a crossing of the phantom divide by using the reconstruction method and investigate the 
behavior of F{R) gravity around a crossing of the phantom divide. Finally, conclusions are 
given in Sec. IV. 



II. RECONSTRUCTION METHOD 



First, we explain the reconstruction method of F{R) gravity proposed in Ref. 33|]. The 
action of F{R) gravity with matter is as follows: 

-F{R) 



S = d Xy/—g 



„ o ~r ^matter 



(2.1) 



where g is the determinant of the metric tensor g^,^ and /^matter is the matter Lagrangian. 
By using proper functions P{(f)) and Q{(f)) of a scalar field 0, the action in Eq. (12. ip can be 
rewritten to 



S= J rf'xv^|^[P(0)i? + Q(0)]+A,attcr| . (2.2) 

The scalar field may be regarded as an auxiliary scalar field because it has no kinetic term. 
From Eq. (12.11) . the equation of motion of (f) is given by 

= ^^^fi+^, (2^3) 

d<p d(f) 

Substituting </> = (f){R) into the action in Eq. (12. 2p yields the expression of F{R) as 

F{R)=Pi(PiR))R + Qi(f)iR)). (2.4) 
From Eq. (12.21) . the field equation of modified gravity is derived as 



-g,, [Pi<P)R + g(0)] - R,.Pi<P) - g,uaP{<p) + V,V.P(0) + «:2Tj-tter) ^ q ^ ^2.5) 



where is the covariant derivative operator associated with g^^, □ = g'^'^'V^Vi, is the 
covariant d'Alembertian for a scalar field, and j'^™'***"'') jg the contribution to the energy- 
momentum tensor form matter. 

We assume the fiat Friedmann- Robert son- Walker (FRW) space-time with the metric, 

ds^ = -dt^ + a^{t)dx^ , (2.6) 
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where a(t) is the scale factor. In this background, the components of (/i, z/) = (0, 0) and 
(/i, i^) = (z,j = I,-- - ,3) in Eq. ([23]) read 



SH^P{(t){t)) - Q{(t){t)) - 6H 



dt 



+ 2k^p = 



(2.7) 



2^lI^M + AH^^M^ +Uh + 6hA P(0(t)) + g(0(t)) + 2k'p = , (2.8) 



dt^ dt 

where H = a/ a is the Hubble parameter with ' = d/dt and p and p are the sum of the 
energy density and pressure of matters with a constant EoS Wi, respectively, with i being 
some component of matters. After eliminating (5(0) from Eqs. fl2.7p and (12. 8p . we get 

d-'pm)) ^dpm) 



2HP{(f){t)) + K' {p + p) = 0. 



(2.9) 



df^ " dt 

The scalar field (f) may be taken as = t if it is redefined properly. By representing a{t) as 



a{t) = aexp {g{t)) 



(2.10) 



in terms of a constant of a and a proper function of g{t) and using H = dg{(f))/ {d(j)), we 
rewrite Eq. (12. 9p to be 

rf2p(0) d~g{<P)dPi^) , 



+ 2- 



-P{<f>) 



where pj is a constant. Moreover, from Eq. (12. 7p . we obtain 



(2.11) 



Q{4>) = -6 



d^(</)) 



P(0) - 6 



dm dP{<P) 



+ 2k^ Y1 Pi«"'^'^"'^ exp [-3 (1 + w,) g{(j))] . 



(2.12) 



We note that if we redefine the auxiliary scalar field by = $(</?) with a proper function 
$ and define P{ip) = P{^{ip)) and Q{ip) = Q{^{(p)), the new action 



S 



d'^x^f—g 



F{R) 



2k^ 



F{R)=P{^)R + Q{^) 



(2.13) 
(2.14) 



is equivalent to the action in Eq. (12.21) because F{R) = F{R). Here, cp is the inverse function 
of $ and we can solve ip with respect to R as ip = 'p{R) = ^~^{(j){R)) by using = 0(i?). 
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Consequently, we have the choices in hke a gauge symmetry and therefore we can identify 
(p with time t, i.e., (p = t, which can be interpreted as a gauge condition corresponding to 



the reparameterization of = (f){(p) [26|. Thus, if we have the relation t = t{R), in principle 
we can obtain the form of F{R) by solving Eq. (12.111) with Eqs. (12. 4p and fl2.12p . 

We also remark that a crossing of the phantom divide cannot be described by a naive 
model of F{R) gravity. To realize the crossing, F{R) needs to be a double-valued function, 
where the cut could correspond to Wcs = — 1. However, the crossing can be performed by 
the extension of F{R) gravity, whose action is given by P{(f))R + Q{(f))- 

III. MODEL 

Next, we examine a model of F{R) with realizing a crossing of the phantom divide by 
using the reconstruction method and investigate the behavior of F{R) gravity around a 
crossing of the phantom divide. 

A. Crossing of the phantom divide 

To illustrate the behavior of F{R) with realizing a crossing of the phantom divide, we 



consider the case in which the Hubble rate H{t) is expressed as 34 1 



ff = Ht + r^l. (3.1) 



35| appears 



where n is a positive constant and ts is the time when the Big Rip singularity 
as will be shown later^. 

In the FRW background, the effective energy density and pressure of the universe are 
given by pcs = SH'^/k'^ and Pefr = — (^2H + 3H^^ /k^, respectively. The effective EoS 



Wcs = Pes/PeS is defined as {g] 



^efr = -l-|^, (3.2) 



which implies that a crossing of the phantom divide occurs when the sign of H changes. 

When t — i> 0, i.e., t <^ tg, FI{t) in Eq. ( 13. ip behaves as H{t) ~ n/t and therefore 
H ~ —n/t^ < 0. In this limit, it follows from Eq. (13. 2p that the effective EoS is given by 



^ Other kinds of finite-time future singularities have been studied in Ref. 36 1 
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FIG. 1: Time evolution of ifJefr for n = 10 and a = 2.0 with t = t/tQ. 



Wes = —1 + 2/ (3ra) > —1. Such behavior is identical with that in the Einstein gravity with 
matter whose EoS is greater than —1. This is the non-phantom phase. On the other hand, 
when t ts, we find H{t) ~ n/ (tg — t) and hence H n/ (t^ — t)^ > 0. We only consider 
the period < t < tg because H should be positive. In this case, the scale factor is given by 
a{t) a {tg — Thus, when t — > t^, a — > cxo, namely, the Big Rip singularity appears. In 
this limit, the effective EoS is given by Wcs = —1 — 2/ {3n) < —1. Such behavior is identical 
with the case in which there is a phantom matter with its EoS being smaller than —1. This 
is the phantom phase. Moreover, from Eq. (13. 2p . we see that the effective EoS Wce becomes 
— 1 when H = 0. Solving Wes = — 1 with respect to t by using Eq. (13. ip . we find that the 
effective EoS crosses the phantom divide at t = tc given by tc = ts/2. As a consequence, in 
case of Eq. (13. ip . a crossing of the phantom divide can occur. We show the time evolution of 
Wcs in Fig. 1 with t = t/tQ, where to is the present time. In all figures, we take n = 10 and 
ts = ato with a = 2.0. In this case, tc = to- From Fig. 1, we see that at the present time, 
a crossing of the phantom divide from the non-phantom phase {Wcs > —1) to the phantom 
one {wcfi < — 1) can be realized. 
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B. Behavior of F{R) gravity around a crossing of the phantom divide 



In what follows, we take (j) = t. From Eqs. (12.101) and (13. ip . H = dg{t)/ (dt) and 
R = q(h + 2hA , we obtain 



g{t) = rain 
a{t) 
R 



t 



ts-t 

a-l)t 



ts-t 
6nt^ 



[{2n -l)ts + 2t\ 



(3.3) 
(3.4) 
(3.5) 



(t. - tf 

where we have taken a = {a — 1)" so that the present value of the scale factor should be 
unity. 

We define X = t/ts and solve Eq. (13.51) with respect to X. If n is much larger than unity, 
we can neglect the second term on the right-hand side of Eq. (13. 5p and therefore obtain the 
approximate solutions 



X{R) 



1 ± Vl - AR-^I^ 



where 



R 



tiR 



(3.6) 



(3.7) 



6n (2n - 1) 

For the lower sign and the upper one in Eq. (13. 6p . X varies asO<X<l/2 and 1/2 < X < 1, 
respectively. 

For simplicity, we consider the case in which there exists a matter with a constant EoS 
w = p/p. In this case, Eqs. (12. lip and (I2.12p are rewritten to 

dR"" 



AR^I- fl- 4^-/2) 

V / WR2 



2R' 



(?> - 10^-^/2j ^ n^/l-m-y^^ ± 2nRVl - 4R-y^P{R) 



+ t;K\i + w)p 



(a-l] 



1 ± Vl -4^-^/2^ 
1 T Vl - 4:R-y'^ 



-3n{l+w) 



(3.8) 



and 



isQiR) 



n 



2n-l ' 2n-l dR 



+ 



3n {2n - 1 



.2 2 - 
-t^K p 



ia-1] 



1 ± VI -4i?-i/2 
1 ^ - 4^-1/2 



-3n(l+ii>) 



(3.9) 



respectively, where ts = tg/ [Qn {2n — 1)]. Here, p corresponds to the present energy density 
of the matter. In particular, we use the present value of cold dark matter with w = for 



p, i.e., p = 0.233pc [H, where pc = 
density and Hq = 2.13h x 10~^^GeV |37j with h = 0.70 |38|, [sgj is the present Hubble 
parameter. From Eqs. fl2.4p and fl3.7p . we have 

We examine F{R) by solving Eqs. fl3.8p - fl3.10p numerically. 

In Figs. 2 and 4, we depict P{R) and t^Q{R) as functions of R. We have used to ~ 
I/Hq. Figs. 2 and 4 show the case of < t < tc{= to), i-e., < X < 1/2 and that of 
^c(= to) < t < tg, i.e., 1/2 < X < 1, respectively. We have numerically solved Eq. (13. 8p 
in the range of R as 16.0001 < R < 18.0. Here, we have taken the initial conditions as 
P{R = 16.0001) = 1.0 and dP{R = 16.0001)/(rf^) = so that around the time when 
a crossing of the phantom divide occurs tc, F{R)/ (2k^) could contain the term R/ (2k^), 
namely, the ordinary Einstein- Hilbert action. We note that at t = tc? i-e., X = 1/2 and 
hence R = 16.0, we cannot solve Eq. ( \3.8\} numerically. We therefore investigate the behavior 
of F{R) gravity around a crossing of the phantom divide. By using Eq. f l3.10p . we show 
the behavior of F{R)/ (2k^) in Figs. 3 and 5, which show the case of < t < tc(= to), i.e., 
< X < 1/2 and that of tc(= to) < t < tg, i.e., 1/2 < X < 1, respectively. For a = 2.0, 
because tc = ^o? < X < 1/2 and 1/2 < X < 1 correspond to the past and the future, 
respectively. Furthermore, we illustrate WesiR) = — 1 + [2/ (3n)] ^1 — 2X(i?)j in Fig. 6. 
The time evolution of R is given in Fig. 7. From Figs. 6 and 7, we see that a crossings of the 
phantom divide can be realized. The results in the all figures are shown by dimensionless 
quantities. 

From Figs. 3 and 5, we see that before a crossing of the phantom divide, F{R) decreases in 
terms of R, while that after the crossing, F{R) increases in terms of R. The latter behavior 



(SnG) = 3.97 X 10-^^GeV^ is the critical energy 



is reasonable because in the Hu-Sawicki model [40|] of F{R) gravity, which passes the solar 
system tests, F{R) increases around the present curvature. For a = 2.0, the time when a 
crossings of the phantom divide is the present time, and it follows from Eqs. (13. 5p and (13. 7p 
that the present value of R is Ro = 16.0. We remark that such behavior is typical for a 
general class of viable F{R) gravities introduced in Ref. 4l| to which the Hu-Sawicki model 



belongs. This class of modified gravities can satisfy the solar system tests and unify inflation 
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FIG. 2: P{R) and i^Q{R) as functions of ^ for < t < tc{= to). We have taken n = 10, a = 2.0 
and p = 0.233pc- 
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FIG. 3: Behavior of F{R)/ (2k^) as a function of i? for < t < tc{= to). Legend is the same as 
Fie. 2. 
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FIG. 4: P{R) and i^Q{R) as functions of R for tc{= to) < t < tg. We have taken n= W, a = 2.0 
and p = 0.233pc- 




16 16.5 17 17.5 18 



R 

FIG. 5: Behavior of F{R)/ as a function of for tc(= to) < * < ^s- Legend is the same as 

Fig. 4. 
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FIG. 6: Behavior of Wcs{R) as a function of R. The left panel and the right one show the case of 
< t < tc{= to), i.e., < X < 1/2 and that of tc{= to) < t < tg, i.e., 1/2 < X < 1, respectively. 
Legend is the same as Fig. 1. 
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FIG. 7: Time evolution of R. Legend is the same as Fig. 1. 
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with the late-time cosmic acceleration. As viable models of F{R) gravity, e.g., the models 



in Refs. 



41 



42 



43 



44 



45 



46 



47| are also proposed. In addition, investigations to solve the 



problem of a curvature sin gula rity in F(R) gravity 



54 



55 



56 



57 



58 



17 



48 



49 



50 



51 



52l | have recently 



59| . Theories without a singularity were also 



been executed in Refs. |53|, 
constructed in Ref. jl" 

Finally, we mention the stability for the obtained solutions of the phantom crossing under 



a quantum correction coming from the conformal anomaly. In Ref. [26|], it has been shown 
that the quantum correction of massless conformally-invariant fields could be small when a 
crossings of the phantom divide occurs and therefore the solutions of the phantom crossing 
could be stable under the quantum correction, although the quantum correction becomes 
important near the Big Rip singularity. 



IV. CONCLUSION 



In the present paper, we have investigated a model of F{R) gravity in which a crossing 
of the phantom divide can be realized. In particular, we have illustrated the behavior of 
F{R) gravity around a crossing of the phantom divide by taking into account the presence 
of cold dark matter. The demonstration in this work can be interpreted as a meaningful 
step to consider a more realistic model of F{R) gravity, which could correctly describe the 
expansion history of the universe. 
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